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We summarize the status of determinations of fundamental parameters of QCD by the ^isESi. 



1. INTRODUCTION 

The natural choice for the fundamental param- 
eters of QCD are the renormalization group in- 
variant (RGI) quark masses, M^i = u, d, and 
the A-parameter. The RGI masses do not depend 
on a renormalization scheme and the scheme de- 
pendence of the A-parameter is exactly given by a 
one-loop coefficient. These quantities are in prin- 
ciple free of perturbative uncertainties. 

Their determination from low energy hadron 
data is a relevant mission for lattice QCD. Before 
entering details, we discuss a few aspects of more 
conventional determinations of a (and A). One 
considers experimental observables Oi depending 
on an overall energy scale q and possibly some ad- 
ditional kincmatical variables denoted by y. The 
observables can be computed in a perturbative 
series which is usually written in terms of the MS 
coupling a m = g^/^ir, 2 (n = q) 

Oi{q,y)=a m ( l i)+A i {y)a^ § { l i) + ... . (1) 

For example Oi may be constructed from jet cross 
sections and y may be related to the details of the 
definition of a jet. 

The renormalization group describes the en- 
ergy dependence of g in a general scheme, 
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here the /3-function has an asymptotic expansion 

m -g 3 {b + g 2 b 1 + ...} , 

60 = ^(11-1^) , 
*i = ^(l02-fJV f ) , (3) 

with higher order coefficients bi, i > 1 that de- 
pend on the scheme. 

We note that - neglecting experimental uncer- 
tainties - c^Mg extracted in this way is obtained 
with a precision given by the terms that are left 
out in eq. (1). In addition to (say) a 3 -terms, 
there are non-perturbative contributions which 
may originate from renormalons, instantons or, 
most importantly, may have an origin that no 
physicist has yet uncovered. Empirically, one ob- 
serves that values of determined at different 
energies and evolved to a common reference point 
using eq. (2) including bi agree rather well with 
each other; the aforementioned uncertainties arc 
apparently not very large. Nevertheless, the stan- 
dard determinations of a are limited in precision 
because of these uncertainties and in particular if 
there was a significant discrepancy between a de- 
termined at different energies and/or processes it 
would be difficult to say whether this was due to 
the terms left out in eq. (1) or was due to terms 
missing in the Standard Model Lagrangian. 

It is an obvious possibility, and at the same 
time a challenge, for lattice QCD to achieve a de- 
termination of a in one non-perturbatively (NP) 
well defined scheme at large /1. There one may use 
the perturbative approximation for (3(g) inserted 
into the exact solution of eq. (2), 

A = M (6 5 2 (M))" 61/(26§) e- 1/(2bo52W) (4) 
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in order to extract the A-parameter. 

1.1. Problem 

In an implementation of this idea one needs to 
satisfy the following criteria. 

• Compute a{p) at energy scales of fi ~ lOGeV 
or higher in order to connect in a controlled man- 
ner to the perturbative regime. 

• Keep jU removed from the lattice cutoff a -1 to 
avoid large discretization effects and to be able 
to extrapolate to the continuum limit. 

• Avoid finite size effects in the MC simulations 
by keeping the box size L large compared to the 
inverse of the mass of the pion, m,. 
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Figure 1. The strategy for a non-perturbative 
computation of short distance parameters. 



These conditions are summarized by 

L » i ~ orikv >y Jl~ Tokv >>a ' 

which means that one must perform a MC- 
computation on an A^ 4 lattice with N = L/a ^> 
70: an impossible task for some time to come. 

1.2. Solution 

The simple solution is to identify the renormal- 
ization scale and the box size L [1], and step up 
the energy scale recursively. Each step then re- 
quires only L/a^> 1, see Fig. 1. 

2. RUNNING OF THE COUPLING 

We choose a NP defined Schrodinger functional 
renormalization scheme (SF) and set all renor- 
malization conditions at zero quark mass (see [2] 
and references therein). An essential property of 
this scheme is the use of Dirichlet boundary con- 
ditions in time, which allow for MC simulations 
at vanishing quark mass. 

2.1. The step scaling function 

One vertical step in the figure is achieved as fol- 
lows: we start from a given value of the coupling, 
u = g 2 (L). When we change the length scale by 
a factor 2, the coupling has a value g 2 (2L) = u' . 
The step scaling function, a, is then defined as 

a(u) = u' . (5) 



It is a discrete /3-function whose knowledge al- 
lows for the recursive construction of the run- 
ning coupling at discrete values of the length 
scale, Uk — .g 2 (2~ fc £ max ) , once a starting value 
itmax = <? 2 (£max) is specified. The step scaling 
function has a perturbative expansion u(u) = 
u + 2b ln(2)u 2 + ... . 

On a lattice with finite spacing, a, it has an 
additional dependence on the resolution a/L. Its 
computation then consists of the following steps: 

1. Choose an (L/a) 4 lattice (1 < x^/a =< L/a). 

2. Tune the bare coupling g such that the renor- 
malized coupling g 2 {L) has the value u and tune 
the bare quark mass too such that the PCAC- 
mass vanishes. 

3. At the same values of go, mo, simulate a 
(2L/a) 4 lattice; compute u' = g 2 (2L). The lat- 
tice step scaling function is S(u, a/L) = v! . 

4. Repeat steps 1.-3. with different resolutions 
L/a and extrapolate ct(u) = lim a /£_^ Y,(u,a/L). 

Note that step 2. takes care of the renormal- 
ization and 3. determines the evolution of the 
renormalized coupling. The continuum extrapo- 
lation 4. is an essential part. In QCD, starting 
with the pure gauge theory, a lot of effort - both 
perturbative, analytical and by MC - has been 
invested in understanding this step properly. We 
do not have the space here to explain it in de- 
tail, but as a result of [4,5,3,2,6] and other work 
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Figure 2. The QCD /3-function in the SF scheme. 



cited in these references, we are convinced that 
no systematic lattice errors are left. 

3. RESULTS 

We proceed directly to discussing the results. 
Numerical data for the step scaling function have 
been obtained for 7 values of u in the range of 
u = 1.0 to u = 3.3 for QCD with Nf = 2 dynam- 
ical flavours [3,2], while for Nf = there are 11 
points in u between u = 0.9 and u = 2.4 [7,5]. 
After taking the continuum limit, the extracted 
<j(u) have a precision of around 0.5% at the small- 
est coupling (where it runs slowly) to 2% at the 
largest one. The numerical values of cr(u) are then 
represented by a smooth interpolating function (a 
polynomial in u). With this function the running 
coupling can be constructed recursively; the re- 
sult is shown in Fig. 3 and Fig. 4. 

With a start value for the /3-function taken 
from perturbation theory (3-loop) at the weak- 
est coupling (a » 0.08), one can also set up a 
recursion for the /3-function itself [2]. As Fig. 2 
shows, its agreement with perturbation theory is 
excellent at weak couplings a < 0.2, while at the 
largest couplings significant deviations from per- 
turbation theory are present for Nf = 2. Indeed 
the difference between non-perturbative points 
and 3-loop can't be described by an effective 4- 
loop term with a reasonable coefficient. At the 
same time the perturbative series just by itself 



does not show signs of its failure at, say, a « 0.4. 
Similar deviations from perturbation theory, at 
even somewhat smaller values of a have been 
found in the NP scale dependence of composite 
operators [8,9]. 

We return to the running couplings shown in 
Fig. 3 and Fig. 4. In the zero flavour case, also the 
region of ^ of around 250 MeV was investigated 
with a specifically adapted strategy [10]. In this 
region, the SF coupling shows the rapid growth 
expected from a strong coupling expansion. 

Initially, the graphs Fig. 3 and Fig. 4 are ob- 
tained for /i in units of /i m in = 1/ imax defined by 
3 2 (£ max ) = Umax. One chooses u max relatively 
large, but within the range covered by the non- 
perturbative computation of o~(u). Since it has 
been verified that perturbative running sets in (in 
the SF scheme) at the larger il that were simu- 
lated, one can use eq. (4) with the perturbative 
/3-function in that region. In this way L max can 
be expressed in units of A -1 : 

Nf = u max = 3.48 

ln(A m L max ) = -0.84(8) (6) 

Nf = 2 u max = 4.61 

ln(A m L max ) = -0.40(7) . (7) 

It remains to relate the artificially defined length 
scale L max to an experimentally measurable low 
energy scale of QCD such as the proton mass or 
the kaon decay constant, Fk- So far we have 
evaluated L max in units of the scale parame- 
ter ro, which has a precise definition in terms 
of the QCD static quark potential [11] and is 
well calculable on the lattice. But it is related 
to experiments only through potential models: 
ro ~ 0.5 fm. For Nf = the continuum limit 
is L max /r = 0.738(16), A^r = 0.60(5) [12]. 



Table 1 

Lambda parameter in units of ro for Nf = 2. 





^max 


= 4.61 


r /a 




A MS r 


5.45(5)(20) 


6.00(8) 


0.610(25) 


6.01 (4) (22) 


6.57(6) 


0.614(24) 


7.01(5)(15) 


7.73(10) 


0.609(16) 
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In the N{ = 2 theory, the situation is illus- 
trated in Table 1 which relies on results for r^/a 
from [13,14]. On the one hand, all the numbers 
are consistent, indicating that lattice spacing ef- 
fects are small, on the other hand the first colum 
shows that a is not yet varied very much. At the 
moment we quote Ai-jLj"o = 0.62(4)(4), where the 
second error generously covers the range of Ta- 
ble 1 as well as numbers obtained with u max = 
3.65 [2] and the first one comes from eq. (7). 

4. DISCUSSION 

The scale dependence of the SF coupling is 
close to perturbative below asF = 0.2 and be- 
comes non-perturbative above o;sf = 0.25. In 
fact a strong coupling expansion suggests that 
this coupling grows exponentially for large L and 
in the Nf = theory it was possible to verify this 
behavior for L close to lfm [10] (Fig. 4). 

In the following table, we compare our results 
for Aj^g-ro to selected phenomenological ones, 
where we set r — 0.5 fm. 

N { : 2 4 5 

[2,7] 0.60(5) 0.62(6) 

[15] 0.74(10) 0.54(8) 
[16] 0.57(8) 

There appears to be an irregular in- 
dependence, but we note that 

1. the errors are not very small, 

2. the 4-flavour A is related to the 5-flavour one 
by perturbation theory [17] which is then required 
to be accurate for \x <C mbeauty- 

4.1. Improvements of lattice results 

The replacement of L max /r by F K x L max 
(computed at small enough light quark masses 
and small a) is the most urgent step and an ef- 
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Figure 3. Running coupling for Nf = 2. 



fort is presently being made. Also the strange 
quark sea has to be included ("2+1") and one 
should estimate the effects of the charm quark. 
Such 2+1 simulations are for example being car- 
ried out by JLQCD and CP-PACS, who have also 
studied the computation of the SF coupling with 
the gauge action they are using [18]. 

I am convinced that lattice results will yield 
the best controlled and most precise A in the long 
run. The reason is simple: I essentially described 
all the sources of systematic errors and the kind 
of assumptions one has to make are minimal. 

There is a large number of other results from 
lattice gauge theory, where a(fj,) is extracted from 
a quantity related to the cutoff: /i ~ a" 1 . Some 
of them cite a very small error [19]. As discussed 
for example in [20], it is, however, not easy to 
estimate the systematic errors in these computa- 
tions, mainly because one cannot separate higher 
order perturbative corrections from discretization 
errors. It thus remains very desirable to carry out 
the program of Fig. 1 with good precision and the 
relevant number of flavors. 

5. QUARK MASSES 

The strategy discussed above is easily gener- 
alized to include the quark masses [5]. A useful 
ingredient is again to choose a mass independent 
renormalization scheme (even for mt,) where (i = 
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Figure 4. Running coupling for Nf = 0. 
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Table 2 

Quark masses determined with full NP renormal- 
ization and continuum limit. We use r = 0.5 fm. 



i 


N { 


input 


Mi/GeV 


ref. 


strange 





m K ,r 


0.137(05) 


[22] 


strange 


2 


m K ,r 


0.137(27) 


[21] 


charm 





m D , r 


1.654(45) 


[23] 


beauty 







6.771(99) 


[24] 



up, down, strange ...) 

m i (l*)=Z m ( f M,g )m% m (g ) (8) 

with a flavour independent rcnormalization factor 
Z m {fia, go). The scale dependence of the running 
masses can then be computed once and for all. 
The RGI masses are defined by 

Mi = lim (26 5 2 ( M ))- d °/ 2bo m( M )i (9) 

fj, — *oo 

and we have 

Mi/Mj = mitii/mitij = m^/mf™ . (10) 

In [5,21] the renormalization problem has been 
solved non-perturbatively in the Nf — 0, 2 the- 
ories: the /z-dependence in the SF scheme was 
computed at the % level. Based on this the re- 
sults in Table 5 were obtained. They may easily 
be converted to rny^(fi) if the desired renormal- 
ization scale n is not too small. 

For the charm and the beauty quark masses 
determinations with Nf > and NP renormaliza- 
tion are still missing. However, in our opinion it 
is even quite early concerning the determinations 
of the light quark masses. Although some Nf- 
dependence of the strange quark mass has been 
reported in the literature, one can presently not 
exclude that this is due to pcrturbative uncertain- 
ties or discretization errors. 

Note that the quark masses computed in the 
quenched approximation were in a similar stage 
in 1996 but very soon afterwards the uncertain- 
ties shrunk by an order of magnitude due to NP 
renormalization and continuum extrapolations! 
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